Annihilation Type Radiative Decays of B Meson in Perturbative QCD Approach 
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With the perturbative QCD approach based on kr factorization, we study the pure annihilation 
type radiative decays B° — > 07 and B° — > J/ipj. We find that the branching ratio of B — > (fry is 
(2.7lQgln'fi) x 10 -11 , which is too small to be measured in the current B factories of BaBar and 



Belle. The branching ratio of B° —> J/ipl is (4.5 



+0.6+0.7\ 



x 10 , which is just at the corner of 



being observable in the B factories. A larger branching ratio BR(B® — > J/ip~f) — 5 x 10~ 6 is also 
predicted. These decay modes will help us testing the standard model and searching for new physics 
signals. 
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PACS numbers: 13.25.Hw, 12.38.Bx 

B meson rare decays are interesting for testing the 
standard model and searching for new physics. However, 
due to our poor knowledge of non-perturbative QCD, 
predictions for many interesting exclusive decays are pol- 
luted by large hadronic uncertainties. The two body ra- 
diative B decays involve simpler hadronic dynamics with 
only one hadron in the final states, so they suffer much 
less pollution than non-leptonic decays. The radiative 
decays such as B — > K*j, p(u>)j thus attract much atten- 
tion [1] . The isospin breaking effects between the charged 
B^ and neutral B° in these modes are mainly due to con- 
tributions from the annihilation type diagrams [1-5]. 

The importance of the annihilation type diagrams can 
also be shown from the pure annihilation radiative B de- 
cays. The color suppressed B° — ► J/ipj and B° — ► <fi-f 
modes are of this kind. The former is tree dominant 
while the latter is a pure penguin flavor changing neutral 
current decay. Despite the fact that they are annihila- 
tion type decays, these decay amplitudes can be factor- 
ized as the B meson to photon transition form factor 
(tI?Tm(1 ~ 75)^1^) times the decay constant (0|s , y M s|O) 
or (J/tplcjfiClO) in the naive factorization approach. 

Recently, the vertex corrections for the four quark 
operators have been performed in the so called QCD 
factorization approach [6], utilizing the light-cone wave 
functions. The branching ratios turn out to be one or- 
der of magnitude different from the naive factorization 
approach [7, 8]. Such a large contribution from next- 
to-leading order corrections implies that the hadronic 
uncertainty in this kind of decays is as large as other 
hadronic annihilation type decays [9]. Recent study of 
soft collincar effective theory [10] also shows that the 
naive factorization contribution is not the only domi- 
nant contribution, which contradicts to QCD factoriza- 
tion claims. More theoretic study is needed before one 
can claim new physics effects in these decays. 

In this paper, we will use an alternative approach- the 
perturbative QCD approach (PQCD) [11] to calculate 



the pure annihilation type decays B a — > ^7 and B° — > 
J/ip"f. Based on k^ factorization, the PQCD approach 
has been proposed and applied to calculate two body 
non-leptonic B decays such as B — > Kir [12], 7T7t [13], pp 
[14], the radiative decays B — > K*j [4], p(ll>)j [5], etc. 
and the results are consistent with experimental data. 
In this approach, the quark transverse momentum ky is 
kept in order to kill the end-point singularity. Because 
of inclusion of transverse momenta, double logarithms 
from the overlap of two types of infrared divergences, 
soft and collinear, are generated in radiative corrections. 
The resummation of these double logarithms leads to a 
Sudakov form factor, which suppresses the long-distance 
contribution. 

For convenience, we work in the light-cone coordinate, 
where the B meson momentum in its rest frame, is 



Pb — (Pb, p b ' Pe 



Mb 
V2 



(1,1, Ox). 



(1) 
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By choosing the coordinate frame where the vector me- 
son moves in the "— " direction and photon in the "+" 
direction, the momenta of final state particles are 

P V = (P+,Py,P v± ) = ^(r 2 ,l,6 ± ) 7 

P-y = (P+,P-,P 7 x) = ^(l-r 2 ,0,0x), (2) 

where r = my /Mb ■ The momentum of the light quark 
in B meson is: 

, _ - Mb 
ki = {k{ , fe x ,k\_ T ) = (— f=-x\,Q,kiT)- (3) 

For the final state vector meson, we set the momentum 
of q(q = s, c) as 

fc 2 = (fc^,fc^,fc 2T ) = (^^X2r 2 ,^^x 2 ,k 2 T)- (4) 

In above functions, x\ and 22 are momentum fractions 
of the quarks. 

In PQCD approach, the decay amplitude is factorized 
into the convolution of the mesons' wave functions, the 
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FIG. 1: Feynman diagrams for B° — > 07 process in PQCD. 



hard scattering kernel and the Wilson coefficients , which 
stand for the soft, hard and harder dynamics respectively. 
With transverse momentum and Sudakov form factor, 
the formalism can be written as: 

/>1 poo 

/ dx x dx 2 I b x db x b 2 db2 ®v(%2, &2)$b(zi, b x ) 
Jo Jo 

x C(t)H(xi,X2,bi,b 2 ,t)eKp[-S(xi,X2,bi,b2,t)}, (5) 

where bi is the conjugate space coordinate of the trans- 
verse momentum kiT, which represents the transverse in- 
terval of the meson, t is chosen as the largest energy 
scale in the hard scattering kernel H in order to suppress 
higher order corrections. The light cone wave functions 
of mesons are not calculable in principle in PQCD, but 
they are universal for all the decay channels. So that 
they can be constrained from the measured other decay 
channels, like decays B — > Ktt [12], B — > irir [13], etc. 

Since the outgoing photon can be only transversely 
polarized, the decay amplitude can be decomposed into 
two parts as: 



A = (sy ■ e*)M s + 



e^ pa e^e^PPP v M p , (6) 



P v ■ P^""^ 



where Py and P 7 are the momenta of vector meson, and 
photon, respectively. e v and e* are the relevant polar- 
ization vectors. The matrix element M s ^ can be cal- 
culated in the PQCD approach. 

In principal, the leading order contributions for B° — > 
V7 (V = J/ip,<f>) decays involve only four-quark op- 
erators plus a photon emitted from any of the quark 
line. The effective weak Hamiltonian is formed by the 12 
four-quark operators and the corresponding QCD cor- 
rected Wilson coefficients [15]. Very recently the elec- 
tromagnetic penguin operator C>7 7 contribution through 
B° — > 77 with one photon connecting to the <j> meson is 
studied in ref.[16]. The branching ratio for B — > (fry is 
found to be 1 x 10 -11 which is larger than the four quark 
operator contribution from QCD factorization approach 
[7, 8]. In PQCD language, the contribution from Oi 1 
is next-to-leading order. Its contribution is still smaller 
than other contributions in PQCD approach which will 



be shown later. The contribution of this kind of operator 
to the B° — ► J/tpj decay is negligibly small. 

The lowest order Feynman diagrams of B — > (fry in 
PQCD are shown in Fig.l. In principle, the photon can 
be emitted from any quark line of the four-quark oper- 
ator. However, the contribution of Fig. 1(c) is canceled 
exactly by that of (d) because of topology symmetry. 
Thus, only the decay amplitudes from Fig. 1(a) and (b) 
are left as: 

M s = —G F eM 3 B rf vVckm / dx x \ db x b x 

VO Jo Jo 

x 4> B ( Xl , h) [(1 - r 2 )C(t a )^ (&i^a)e- SB(ta) 

-(l-r 2 )C(i b )K (foii?a)e- SB(t6) l, (7) 



M p = -^=G F eM 3 B rf v V C KM / dx x 



db x b x 



1 

v3 jo Jo 

x (j> B {x x ,b x ) [c(t a )K {b x A a )e- SB ^ 

+ C(t b )K (b x B a )e- s ^ 



(8) 



with 



Al = (l+x x - r 2 )M 2 Bl B 2 a =x x (l- r 2 )M 2 B , 



t a = max(Aj, 1/bi), t b = max(B a , l/b x ). 



(9) 



Kq(x) is the modified Bessel function which results from 
Fourier transformation of the quark propagator. e~ SB ^ 
is Sudakov form factor, fy is vector meson decay constant 
and Vckm denotes CKM matrix elements. In above cal- 
culations, the B meson is treated as a heavy-light system, 
whose wave function is defined as: 



(,Pb + M B )j 5 (f) B {x x ,b 1 ), 



(10) 



and the expression of distribution amplitude (f> B is shown 
in rcf. [12, 13] with parameter u>b — 0.4GeV, which is 
normalized as 



dx X (f> B (x x , b x = 0) 



Sb 

2V6 ' 



(11) 



where fs is the decay constant of the B meson. Since 
we do not need the <f> or J /if} wave functions in the calcu- 
lations of Feynman diagrams Fig. 1(a) and (b), we need 
not show them here. 

For B° — > </>7 decay, only penguin operators can con- 
tribute, and the CKM matrix elements are Vckm = 
V t * b Vtd- The combination of the Wilson coefficients in 
eq.(7, 8) are: 



- \c -±c w . 



(12) 
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TABLE I: Summary of input parameters [19] 



CKM parameters and QCD constant 
A A p fj A^T 4) r B o 

0.2196 0.819 0.20 0.33 250MeV 1.54ps 
Meson decay constants 
fs fs s U f.j/4, 

216McV 236MeV 254MeV 405 ± 14MeV 



Meson masses 
M B M 



M 



80.41GeV 5.28GeV 1.02GeV 3.10GeV 



The Wilson coefficient of dominant QCD penguin oper- 
ator C3 cancels much with G\/3, and C5 cancels with 
Ce/3. This is a result of color suppression in this decay, 
since at least three gluons are needed for a vector <fi me- 
son produced from penguin diagram [17]. For B — ► J/ipj, 
both tree and penguin operators give contribution, and 
corresponding conclusions of the Wilson coefficients are: 



°V>7 



Ci + \c 2 , 

c 3 + l -c A + c 5 + l -c, + c 7 

— C& + Cg + gClO- 



(13) 



Again, C\ cancels much with C2/3 which is also a result 
of color suppressed tree contribution. The Vckm in tree 
(penguin) operators is V* b V cd (V* b V td )- 

With the amplitudes M s and M p defined in Eq.(6), 
the decay width of B° — » Vj is given by 



\M S \* + \M P \* 
8ttM b { 



r 2 )- 



(14) 



In our numerical calculations, the input parameters are 
summarized in Table I, where A, A, p and r\ are CKM 
parameters in Wolfenstein paramctrization [18], and p = 
p(l — 5A 2 ), fj = 77(1 — 5A 2 ). Their values can be found 
in Review of Particle Properties [19]. 

At the leading order, the main uncertainty for decay 
branching ratios comes from the B meson wave function. 
But it is constrained by the measured exclusive hadronic 
decays, like B — > Ktt [12], B — > 7T7t [13] with parameter 
ijj B from 0.3GeV to 0.5GeV. On the other hand, there 
should be large uncertainty since we work only at leading 
order in a s for the hard part and also for the Wilson coef- 
ficients. The missing next-to-leading order correction is a 
very important uncertainty for rare decays. To estimate 
it, we consider the hard scale t at a range of 

max(0.75yl a , ^-) < t a < max(1.25A a , ^-), (15) 

max(0.75 J B 6 , -!-) < t b < max(1.25B 6 , -!-). (16) 
bi bi 



With the above major uncertainties from B meson wave 
function parameter ajg and the different scale t, respec- 
tively, we give out the branching ratio of B° — ► (fry: 



BR(B° -> 07) = (2.7 



+0.3+1. 2x lf) - 
0.6-0.6J x iu 



(17) 



There are many other uncertainties in our calculation 
such as decay constants and CKM matrix elements. How- 
ever, the uncertainty induced by above factors is not 
more than 10% [20]. With such a small branching ratio 
(17), this decay is too rare to be measured at the running 
B factories or even in future LHC-b experiment. Our 
result is still much smaller than the recent upper limit 
8.5 x 10~ 7 from experiments [21]. If some new physics 
particles enhance this ratio through tree or loop effects 
[7], it may be measurable in near future experiments. 

After similar calculations, we also give the branching 
ratio of B° — > J/ipj decay: 



BR(B° -» J/V>7) = (4.5 



+0.6+0.7N ln - 
0.5-0.6) x 1U 



(18) 



Because this process is tree diagram dominated, it is 
much larger than that of B° — ► (fry. This result is still 
not big enough to be measured at the B factories but 
it is already around the corner of B factories capability. 
Currently, there is only upper limit, which is 1.6 x 10~ 6 
at 90% confidence level [22]. 

From our calculation, we find that most of the contri- 
bution comes from Fig. 1(b), which indicates that photon 
emission from the light quark of the B meson is easier 
than that from the heavy quark. The reason is that the 
heavy quark is much difficult to become off-shell than the 
light quark, while a nearly on-shell quark rarely emits 
photons. 

Within QCD factorization approach, the branching ra- 
tios of these two decays have been calculated in Ref. [7, 8] . 
The results are given as : 



BR(B° J/^7) = 7.65 x 10~ 9 ; 
BR(B° 7 ) = 3.6 x 10~ 12 . 



(19) 
(20) 



Compared with their results, we find that our results are 
one order of magnitude larger. In QCD factorization ap- 
proach, some next-to-leading order contribution has been 
added in addition to the naive factorization contribu- 
tion. It shows that these presumably order a s corrections 
change the branching ratios by one order of magnitude. 
In fact these two decays are non-factorizable contribu- 
tion dominant in naive factorization approach. Usually 
QCD factorization approach calculation for this kind of 
decays (such as B° -» L»°7r° [23]) is not stable. They 
receive dominant contribution from non-factorizable di- 
agrams, sometimes with an endpoint singularity in the 
collinear factorization. In fact, our numerical results are 
also much larger than the naive factorization approach, 
which is mainly due to a smaller energy scale (^/ Aqc Dm-b 
other than mi,) used in the calculation of the Wilson co- 
efficients. This uncertainty comes from part of the next- 
to-leading order effect in PQCD approach. 
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The B s — > J/^>7 decay is the same as i3° — > J/ijj'y 
decay within SU(3) symmetry. It is also a pure an- 
nihilation type decay. The decay formulas are exactly 
the same except replacing the corresponding CKM fac- 
tor V c d, V t d by V cs , V ts , respectively. Taking into account 
the SU(3) breaking effect, we used uo Bs = 0.5GcV [24], 
fs s — 236MeV, we can get 

BR{B S -» J/tpj) ~ 5 x 1CT 6 , (21) 

which should be easy to be measured in the future LHCb 
experiment. 

In a summary, we calculate the branching ratios of 
pure annihilation type radiative decays B° — > <fij and 
B° — ► J/^7 within the standard model in PQCD ap- 
proach. We find the branching ratio of B° — > 07 is at 
the order of 10 . This small branching can not be de- 
tected in the running B factories of BaBar and KEK, 
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